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Abstract. We present an extension of the Monte Carlo Event Generator WHIZARD which includes 
chargino production at the ILC at NLO. We include photons using both a fixed order and a 
resummation approach. While the fixed order approach suffers from negative event weights, the 
resummation method solves this problem and automatically includes leading higher order correc- 
tions. We present results for cross sections and event generation for both methods and evaluate 
the systematic errors due to soft and coUinear approximations. In the resummation approach, the 
residual uncertainty can be brought down to the per-mil level. 
This is an updated version of the results presented in [Tl[2]. 

PACS. 12.15.Lk Electroweak radiative corrections - 13.66.Hk Production of non-standard model 
particles in ee-l- interactions - 14.80.Ly Supersymmetric partners of known particles 



1 Introduction 

In many GUT models, the masses of charginos tend to 
be near the lower edge of the superpartncr spectrum, 
and they can be pair-produced at a first-phase ILC 
with cm. energy of 500 GeV. The precise measure- 
ment of their parameters (masses, mixings, and cou- 
plings) is a key for uncovering the fundamental prop- 
erties of the MSSM [3]. Regarding the experimental 
precision at the ILC, off-shell kinematics for the sig- 
nal process, and the reducible and irreducible back- 
grounds [4] need to be included as well as NLO cor- 
rections for chargino production at the ILC which are 
in the percent regime. We here present the inclusion 
of the latter. 



2 Chargino production at LO and NLO 

Fixed order approach 

The total fixed-order NLO cross section is given by 

CTtot(s,mg)= (1) 
CTBorn(s) -|- cr^+s{s, AE^, ml) + 0-2^3(5, AE^.m^), 

where s is the cm energy, me the electron mass, and 
A Ery the soft photon energy cut dividing the photon 
phase space. The 'virtual' contribution is the inter- 
ference of the one-loop corrections ^ with the Born 
term. The collinear and infrared singularities are regu- 
lated by me and the photon mass A, respectively. The 
dependence on A is eliminated by adding the soft real 
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photon contribution cTg — /soft o'Born(s) with a uni- 
versal soft factor /soft(— [S]- We break the 'hard' 
contribution <J2^3{s, AE^, m^), i.e., the real-radiation 
process e~e^ XiXjl^ into a collinear and a non- 
coUinear part, separated at a photon acollinearity an- 
gle A9^ relative to the incoming electron or positron. 
The collinear part is approximated by convoluting the 
Born cross section with a structure function 
Ux-Ae^,^) H: 

aY,,c{AE,A9,s) = 

/ dx^d^2Mxi)\Mbf{x^,s), (2) 

Jae,o 

where Xi denotes the momentum fraction of the re- 
spective incoming beam after photon radiation. The 
non-coUinear part is generated explicitely using exact 
three particle final state kinematics. 

The total fixed order cross section is implemented 
in the multi-purpose event generator WHIZARD ^i9j us- 
ing a 'user-defined' structure function and an effective 
matrix element 

IXeffP - (3) 
(l + MAE^, A))|MbornP + 2i?e(A^bornA^:irt(A)) 

which contains the Born part, the soft-photon factor 
and the Born-1 loop interference term. In the soft- 
photon region this approach runs into the problem of 
negative event weights ^0\ : for some values of 6, the 
2^2 part of the NLO-corrected squared matrix ele- 
ment is positive definite by itself only if AE^ is suffi- 
ciently large, cf Fig.[T] To still obtain unweighted event 
samples, an ad-hoc approach is to simply drop events 
with negative events before proceeding further. 
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Fig. 1. -dependence of effective squared matrix element 
(-^s — 1 TeV). Upper figure: fixed order effective matrix 
element; lower figure: effective matrix element with the 
one-photon ISR part subtracted. Blue/ solid line: Born 
term; red/ dashed: including virtual and soft contribu- 
tions for AE-y = 10 GeV; magenta/ dotted: same with 
AE^ = 0.5 GeV. AO-. = 1°. 
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Fig. 2. Total and relative cross section as a function of ^/s. 
Upper figure: Born (red, "LO"), fixed order (blue, "NLO") 
and fully resummed (green, "NLO") total cross section, 
lower figure: relative fixed order (blue, solid) and fully re- 
summed (green, dashed) higher corrections with respect to 
Born result 



Resummation approach 

Negative event weights can be avoided by resumming 
higher-order initial radiation using an exponentiated 
structure function /isr [11,12J. In order to avoid dou- 
ble-counting in the combination of the ISR-resummed 
LO result with the additional NLO contributions [5], 
we have subtracted from the effective squared matrix 
element the soft and virtual photonic contributions 
that have already been accounted for in CTs+v This 
defines 

IMl'il^ = |A^cffp-2/,oft.ISR|A^BomP (4) 

which is positive for even low AEj cuts for all values 
of 6 (cf Fig. [T|). Convoluting this with the resummed 
ISR structure function for each incoming beam, we ob- 
tain a modified 2-^2 part of the total cross section 
which contains all NLO contributions and in addition 
includes higher order soft and coUinear photonic cor- 
rections to the Born/one-loop interference. This differs 
from the standard treatment in the literature (cf eg. 
[5j) where higher order photon contributions are com- 
bined with the Born term only ("Born-|-"). 

The complete result also contains the hard non- 
collinear 2 — > 3 part convoluted with the ISR structure 
function: 

+ [ dx,dr3f^sR(^i)fmRi^2)\M^^'\^ 
J A(E.e) 

(5) 



The resummation approach eliminates the problem of 
negative weights (cf. Fig. [1]) such that unweighting of 
generated events and realistic simulation at NLO are 
now possible in all regions of phase-space. 

3 Results 

Total cross section and relative corrections 

Figure m shows the cm. energy dependence of the total 
LO and NLO cross section for chargino production for 
the mSugra point SPSla' 3 and the relative correc- 
tions with respect to the Born result. The corrections 
are mostly in the percent regime and can reach 20% 
in the threshold region. 

Cutoff dependencies 

Fig. [3] compares the AE^ dependence of the numer- 
ical results from a semianalytic fixed-order calcula- 
tion with the Monte-Carlo integration in the fixed- 
order and in the resummation schemes. The fixed- 
order Monte-Carlo result agrees with the semianalytic 
result as long as the cutoff is greater than a few GeV 
but departs from it for smaller cutoff values because 
here, in some parts of phase space, |A^effP < is 
set to zero. The semianalytic fixed-order result is not 
exactly cutoff-independent, but exhibits a slight rise 
of the calculated cross section with increasing cut- 
off due to the breakdown of the soft photon approx- 
imation. For AE^ = 1 GeV (10 GeV) the shift is 
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Fig. 3. Total cross section dependence on AE-y-. 'sa' 
(red, dotted) = fixed-order semianalytic result; 'fix' (green, 
dashed) = fixed-order Monte-Carlo result; 'res' (blue, long- 
dashed) = ISR-resummed Monte-Carlo result; (magenta, 
dash-dotted) — same but resummation applied only to the 
2-^2 part. A6-, = 1°. LO: Born cross section. 

about 2permil (5permil) of the total cross section. 
The fully resummed resuh shows an increase of about 
5 permil of the total cross section with respect to the 
fixed-order result which stays roughly constant until 
AE~^ > 10 GeV. This is due to higher-order photon 
radiation. 

For the dependence on the coUinear cutoff A9^, 
the main higher-order effect is associated with photon 
emission angles below 0.1°. For 6^ > 10°, the collinear 
approximation breaks down. 



Ai^^^ leads to corrections in the percent regime. Com- 
paring "Born-|-" and "res-|-", we see that, for y/s > 
SOOGeV, the additional convolution of the interference 
term with /isr is equally in the percent regime and ad- 
ditionally changes the sign of the higher order correc- 
tions. Therefore, all higher order contributions which 
are contained in "res-|-" are significant and have to be 
taken into account. For more details, cf. [T51ll4j . 



4 Conclusions 

We have implemented NLO corrections into the event 
generator WHIZARD for chargino pair-production at the 
ILC with several approaches for the inclusion of pho- 
ton radiation. A careful analysis of the dependence 
on the cuts AEj, AO reveals uncertainties related to 
higher-order radiation and breakdown of the soft or 
collinear approximations. To carefully choose the re- 
summation method and cutoffs will be critical for a 
truly precise analysis of real ILC data. The version of 
the program resumming photons allows to get rid of 
negative event weights, accounts for all yet known 
higher-order effects, allows for cutoffs small enough 
that soft- and coUinear-approximation artefacts are 
negligible, and explicitly generates photons where they 
can be resolved experimentally. Corrections for the de- 
cays of charginos [15ill6j and non-factorizing correc- 
tions are in the line of future work. Equally, a general 
interface to FormCalc |17| in principle allows for the 
inclusion of any lepton-coUider NLO production pro- 
cess in WHIZARD. 



Event distributions 

In Fig.|4]we show the binned distribution of the chargino 
production angle obtained using a sample of unweighted 
events. It demonstrates that NLO corrections to the 
angular distribution are statistically significant and 
cannot be accounted for by a constant K factor. 

Figure [5] shows the difference between the angu- 
lar event distribution resulting from the fixed order 
method and the resummation method. As a reference, 
one standard deviation from the Born event distri- 
bution is also given. It is clearly visible that, in the 
central-to-forward region, even higher order contribu- 
tions are statistically significant. 



dependence of higher order contributions 

Figure [S] shows the magnitude of second and higher 
order photonic effects in different schemes. We com- 
pare the fully resummed version ( "res-|-" ) to the stan- 
dard way of adding higher order photonic corrections 
in the literature ("Born-f"), which consists in only 
convoluting the Born contribution with /isr, and a 
not completely resummed version ( "res" ) where higher 
order photon radiation is taken into account in the 
2 — !• 2 kinematic regions only. Comparison of "res" 
and "res-|-" shows that the additional convolution of 
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Fig. 4. Polar scattering angle distribution for an integrated luminosity of 1 ab~^ at ^/s — 1 TeV. Left: total number of 
events per bin; right: difference w.r.t. the Born distribution. LO (black, dotted) = Born cross section without ISR; fix 
(green, dashed) — fixed-order approach; res (blue, full) — resummation approach. Cutoffs: AE-, = 3 GeV and A9.y = 1°. 
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Fig. 5. Polar scattering angle dependence of difference between events resulting from completely resummed and fixed 
order method: Nrea — Nfi^. Standard deviation from Born (red, solid) is shown 
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Fig. 6. Relative higher-order effects for different methods: (magenta, long dash dotted) = CTres, (blue/ cyan and dash- 
dotted/ dashed) — o"rcs+, and (red, solid) — crBorn+ vs (JBorn. 
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